Topological Entropy of Quantum Hall States in Rotating Bose Gases 
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Through exact numerical diagonalization, the von Neumann entropy is calculated for the Laughlin 
and Pfaffian quantum Hall states in rotating interacting Bose gases at zero temperature in the lowest 
Landau level limit. The particles comprising the states are indistinguishable, so the required spatial 
bipartitioning is effected by tracing over a subset of single-particle orbitals. The topological entropy 
is then extracted through a finite-size scaling analysis. The results for the Laughlin and the Pfaffian 
states agree with the expected values of In y/2 and In \fi, respectively. 
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Quantum Hall states are characterized by topological 
order, in that they can be described by a set of quan- 
tum numbers that are robust against local perturba- 
tions. Examples of such topological invariants include 
the ground-state degeneracy [1], the Chern number [2], 
and the braiding statistics of the quasiparticle excita- 
tions [3] . Interest in the topological underpinnings of the 
quantum Hall effect has surged recently due to the pos- 
sibility of using these states for performing intrinsically 
fault-tolerant quantum computation [3-5]. The poten- 
tial usefulness of certain fractional quantum Hall (FQH) 
states for quantum computation stems from the expec- 
tation that their quasiparticle excitations possess non- 
Abelian fractional statistics, which can be braided in or- 
der to perform topologically-protected logical operations. 
Of the few FQH states that are thought to possess non- 
Abelian excitations (though these alone would be insuf- 
ficient to perform universal quantum computation), the 
most experimentally accessible is the so-called Pfaffian 
state that occurs at a filling factor of v-p = 5/2 in elec- 
tronic FQH systems. For small particle numbers, calcu- 
lation of the overlap between the Pfaffian wavefunction 
and the exact ground state is very good [6-9], though 
some doubt has been recently cast over the validity of 
the Pfaffian wavefunction description [10, 11]. 

FQH states have also been predicted to occur in ro- 
tating neutral Bose gases confined in harmonic traps 
due to the formal equivalence between the Hamiltoni- 
ans describing these systems and two-dimensional (2D) 
electron gases in transverse magnetic fields [12-15]. No- 
tably, the bosonic Laughlin state occurs at a filling factor 
Vb = 1/2 while for the bosonic Pfaffian the filling factor 
is Vb — 1. Here the filling factor is defined as the ra- 
tio between number of particles and vortices. Although 
such states have not yet been experimentally observed 
due to the difficulty in achieving the required high ro- 
tation rates [16], bosonic FQH states have a distinct ad- 
vantage for topological quantum computing in that quasi- 
particle excitations could potentially be simple to excite 
and control [17]. Previous exact diagonalization calcula- 
tions yield a strong overlap between the bosonic Pfaffian 
wavefunction and the exact ground state, though this 
decreases with increasing particle number [18, 19]. It is 



therefore important to calculate global properties such as 
topological quantum numbers in order to provide further 
evidence for the Pfaffian description of the ground state. 

In this work, we focus on one such quantity called the 
topological entanglement entropy 7 [20-22]. Statistical 
mechanics defines the classical entropy as proportional 
to the logarithm of a state's multiplicity. In a similar 
manner, 7 is defined as 7 = hxD, where D > 1 is the 
total quantum dimension of the topological phase [23]. 
A general quantum Hall state in the Laughlin sequence 
with v = 1/p has D = yf>, while the bosonic Pfaffian 
state has D — %/4. (Of course, the v = 1/4 Laughlin 
state is unambiguously distinguished from the Pfaffian 
state by the filling factor) . The topological entropy thus 
provides a powerful tool for classifying different quantum 
Hall states, as long as it can readily be calculated. 

Fortunately, a connection between the von Neumann 
entropy S and topological order has recently been demon- 
strated, from which 7 can be extracted in principle. Sup- 
pose a topologically ordered state is separated into two 
partitions A and B by a circle of radius R. The von Neu- 
mann entropy is defined as S — Sa = — Tr (ftA 1e-Pa)j 
where pa is the reduced density matrix obtained after 
tracing over region B. In large systems, under conditions 
in which B A or vice versa, the von Neumann entropy 
will scale with the length 2-7ri? of the boundary delimiting 
both sections as 



S = a(2nR) - j + 0(1/R), 



(1) 



where a is a nonuniversal coefficient. Cutting the sys- 
tem into multiple subsections and judiciously combining 
the resulting von Neumann entropies, the terms propor- 
tional to the boundary length can be cancelled, leaving 
only the universal topological entropy 7 [20, 21]. In the 
present work, we use exact calculations for small number 
of particles and a finite-size scaling analysis to obtain 7 
for both the bosonic Laughlin and Pfaffian states. 

We consider a zero-temperature gas of bosons confined 
in a cylindrically symmetric harmonic trap that is rapidly 
rotated around the z-axis with a frequency f2. The ax- 
ial trapping frequency uj z is assumed to be much larger 
than that along the radial direction u, so that the gas 
can by considered quasi-2D. This corresponds to the disk 
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geometry that has been used in previous studies of elec- 
tronic FQH systems [24]. Particles interact via a stan- 
dard delta-function pseudopotential whose strength is 
controlled by a 2D coupling constant g = ^/8irhLj£ 2 a/£ z . 
Variables £ — yJh/Mu) and £ z = \J%/Muj z are the char- 
acteristic oscillator lengths along the radial and axial di- 
rections, respectively, and a is the three-dimensional scat- 
tering length. Expressing all lengths in units of I, ener- 
gies in terms of Tilo 1 and frequencies in terms of u, the 
effective 2D Hamiltonian in the frame co-rotating with 
the atoms at frequency is 



H = 



N 



1. 



1 



N 



i<j 



where g = g/£ 2 is the dimensionless coupling constant, 
N is the number of bosons of mass M, L = ^ m, is the 
^-projection of the total angular momentum (which is a 
conserved quantity in this axisymmetric potential), and 
r = (p, 4>) is a particle's position in polar coordinates. 
For the remainder of this work, we set g = 1 . 

Obtaining the topological entanglement entropy con- 
sists of three main steps: evaluation of the Hamil- 
tonian, its diagonalization, and the calculation of the 
ground state's von Neumann entropy. The Hamiltonian 
is expressed in a lowest Landau level (LLL) approxima- 
tion [25] Fock basis of the form 



im>,m,...,^) = n 



(bl) 



■|o>, 



where ¥ m creates a boson with m units of angular mo- 
mentum, and the TV's are the occupation numbers. The 
Hilbert space size is determined by the number of unique 
ways to distribute L units of angular momentum among 
N particles. This number grows very rapidly which limits 
our study to N < 10 for the Laughlin state and N < 13 
for the Pfaffian state. The Hamiltonian is evaluated in 
this basis by writing it in a second quantized form us- 
ing the bosonic field operators ip(r) = ^2 m b m & m (r) ex- 
panded in terms of 2D harmonic oscillator orbitals in 
the lowest Landau level approximation (i.e. zero prin- 
ciple quantum number), $ m (r) = -^p" 1 ^^ / 2 e lm<t> . 
Using these field operators, the Hamiltonian reduces to 



rk b{D ljk i 



ijkl 



where e, = Tiw[N + L(l — ft)] and 

IT 



V 



ijkl 



2 1 1 x l r J v UT. 



H+j,k+l- 



In this Fock basis, the first term in H is equivalent to the 
identity matrix and simply represents an energy offset. 
The problem is then reduced to finding the eigenstates 



of the interaction matrix H lni . These are obtained by 
exact diagonalization using a Lanczos algorithm. 

A particular quantum Hall state is selected by speci- 
fying the total angular momentum according to the rela- 
tionships between L and particle number N established 
in Rcf. 14. The Laughlin state is selected by requiring 
that L = N(N—1) while for the Pfaffian state the relation 
is L = N(N - 2)/2 and L = (N - 1) 2 /2, for even and odd 
./V respectively. As an example, consider the 2-particle 
Laughlin state. It occurs when L = 2(2 — 1) = 2 which 
corresponds to a Hilbert space spanned by the following 
two states: 1 101) = S^J|0) and |020> = -^b\b\\0). Once 
that L is fixed, the Hamiltonian can be constructed and 
diagonalized. This produces a ground state wavefunc- 
tion described in terms of the Fock basis states. In our 
example, the (Laughlin) state vector obtained is 

|* L > = -^(|101>-|020». 

It is straightforward to verify that this corresponds ex- 
actly to the theoretical Laughlin wavefunction 



N 



*L = l[(z i -z j fe-Z i 



M 2 /2<? 2 



i<3 



after expanding in terms of the bosonic field operators 
1 



* L (ri,r 2 , ...,t n ) 



(0|V(ri)V(r 2 )...V(riv)|*L> 



i<j> 



and using z = pe 

Once the ground state wavefunction has been obtained, 
the von Neumann entropy must be calculated. We adopt 
the method of Refs. [26-28] that uses orbital partitioning 
(see Ref. [29] for particle partitioning). This consists of 
separating the state into regions delineated by a specific 
number n (= m max + l) of single particle orbitals, so that 
region A would correspond to m < m max and region 
B to to > TO max . For the two-particle Laughlin state 
example considered above, the reduced density matrix 
PA for n = 1 (meaning only the to = orbital is in 
partition A) is 

Pa = Tr B (|* L )(* L |) 

= | (|l><lj<0|0><l|l> + |0><0|<2|2><0|0» 

= ^(|i)(i| + |o)(o|) 

which has a two- fold degenerate eigenvalue of 1/2. The 
von Neumann entropy is thus pa = —2 x | ln(|) w 0.69. 

Once S is obtained for n — 1, the process is repeated 
numerous times for increasing n, and in principle Eq. (1) 
can be used to obtain the topological entropy 7 from 
a plot of S versus R. What remains to be specified is 
the partition boundary size. The single-particle orbital 
density is ring-shaped and centered at the origin of the 
trap, with radius (p) m = J drp\§ m (p, <f>)\ 2 ~ \Jm + 1 for 
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FIG. 1: Dependence of the von Neumann bipartite entropy S 
with system size R = (p) n _ 1 - The curves S(R) are linear, as 
predicted by Eq. (1), until finite size effects begin to dominate 
for the larger partition sizes. 

large m. We therefore consider the boundary between 
the two regions A and B for a particular choice of n to 
correspond to a circle of size 2n (p) n _ v Since the factor 
of 2tt can simply be combined with a in Eq. (1), we define 
the boundary size when the first n orbitals are kept in 
partition A as 



R = (P)r 



(2) 



A plot of S(R) for N = 5 to 10 is shown in Fig. 1 
for the bosonic Laughlin state. An initial linear increase 
of S(R) is observed, as expected from Eq. (1); however, 
finite-size effects bring S to zero for larger values of R. 
This is because in small- A systems most of the particles 
occupy low-m orbitals whose amplitudes are largest in 
the vicinity of the trap center; the ground-state wave- 
function has a negligeable overlap with higher angular 
momentum orbitals. 

To perform a proper finite-size scaling analysis of the 
small- iV data, we follow the procedure introduced in 
Ref. 26. The value of S(N —> oo) is estimated by plotting 
S as a function of 1 /A for region A containing different 
numbers of orbitals n, as shown in Fig. 2. Our calcula- 
tions were restricted to the range 1 < n < 5 because we 
have too little data for n > 6. To obtain the A — > oo 
values a linear regression was applied. A slight positive 
curvature in the data for the largest values of A might in- 
dicate the emergence of asymptotic values, but the trend 
was not clear enough to enable a more sophisticated anal- 
ysis. The results are plotted in the inset of Fig. 3 and an- 
other linear regression is made to provide the y-intercept. 
We find 7^ = 0.30 ±0.02, which is slightly lower than the 
expected result of 0.35 for the bosonic Laughlin state at 
filling v = 1/2; this might reflect the rather naive linear 
analysis of the finite-size scaling. 

The same procedure was then repeated for the Pfaffian 



FIG. 2: Finite size scaling of the Laughlin state's von Neu- 
mann entropy S with number of particles, for five different 
partition sizes. Linear regressions provide the y-intercepts 
which yield the values of S(N — » oo) shown in the inset of 
Fig. 3. 




FIG. 3: Obtaining the Laughlin topological entropy from 
Eq. (1) and data shown in Fig. 1, after removing data points 
clearly affected by finited size effects. A linear fit on the com- 
bined data produces a y-intercept of —0.57 ± 0.02, leading to 
jl = +0.57±0.02. In the inset, values of S(N — » oo) obtained 
in Fig. 2 are used, increasing the accuracy of the result. A 
weighted linear regression yields 7l = 0.30 ± 0.02 (error bars 
on individual data points are smaller than the data points) . 

state with N = 5 to 13. Unfortunately, a larger amount 
of scatter was present in the Pfaffian version of Fig. 2, 
preventing a clean extrapolation of the S(N) data to the 
large-A limit. Thus, instead of plotting S(N — ► oo) as 
a function of R to obtain 7p, a linear regression is per- 
formed directly on the collection of data points shown in 
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FIG. 4: The topological entropy for the Pfaffian state is ob- 
tained by using the combined data of five to thirteen particles. 
Data points that were obviously influenced by finite size ef- 
fects (for which S(R) decreased with R, as in Fig. 1), were 
removed. A linear regression yields 7p = 0.7 ± 0.1. 

Fig. 4, where n = 1 to 6. This approach produces a topo- 
logical entanglement entropy of jp = 0.7±0.1 consistent 
with the predicted value of ln(-\/4) = 0.69. 

The close correspondence between the Pfaffian numer- 
ical and theoretical results are quite surprinsing consid- 
ering that corrections for finite-size effects have not been 



made. Fitting the aggragate Laughlin results in the same 
manner gives = 0.57 ±0.02 (see Fig. 3) instead of the 
previously obtained jl — 0.30±0.02, suggesting that the 
obtained 7p's accuracy is coincidental. The error ranges 
quoted above strictly reflects the scatter in the numeri- 
cal data, and underestimates the actual uncertainty by 
neglecting systematic errors. In particular, each data 
point is treated as equally valid, whereas the 0(1/ R) 
term in Eq. (1) clearly favors large-i? results. Ideally, 
only the largest-n data would have been kept; however, 
this would not have provided enough data with which to 
extract values of 7. In light of the discrepancy between 
both Laughlin results, a more reasonable estimate of the 
Pfaffian topological entropy is -fp = 0.7 ± 0.3. 

In conclusion, we have calculated the topological en- 
tanglement entropy for both the bosonic Laughlin and 
Pfaffian states for a rotating Bose gas. For the Laugh- 
lin state, we obtain a result of 0.30 ± 0.02 which is al- 
most consistent with the expected value of In \[2 = 0.35. 
For the Pfaffian state we obtain 0.7 ± 0.3. This value of 
the topological entropy is consistent with the expected 
value of ln\/4 = 0.69, though the large amount of scat- 
ter present in Fig. 4 prevents a completely unambiguous 
identification of the state by the Pfaffian wavefunction. 
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